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ABSTRACT 

» ■ 

W e^u se  analytic  models,  valid  for  neutron  scattering,  and 
numerical  calulcations  to  examine  a  method  proposed  by  Drolshagen 
and  Heller  for  computing  atom  diffraction  from  surfaces.  Our 
main  concern  is  Che  manner  in  which  the  use  of  narrow  packets, 
causing  a  partial  coherence  of  the  incident  beam,  might  affect 
the  diffracted  intensity.  We  find  that  Che  demands  of 
representing  beam  coherence  correctly  conflict  with  the 
requirements  for  the  validity  of  the  propagation  scheme.  We 
suggest  some  improvements  that  should  increase  the  reliability  of 
the  method  without  affecting  greatly  its  remarkable  computational 
and  conceptual  advantages. 


I.  INTRODUCTION 

In  several  recent  papers^  Drolsh'agen  and  Heller  have 
presented  a  reformulation  of  the  quantum  theory  of  atom 
scattering  by  surfaces  which  seems  to  have  great  numerical  and 
conceptual  advantages  over  the  traditional  basis  set  expansion 
method.  The  key  idea,  taken  over  from  Heller's  prior  work.  is 
to  write  the  wave  function  of  the  incident  particle  as  a  sum  of 
Gaussian  wave  packets,  to  propagate  them  i n-d e pend en 1 1  y  by  using 
the  time  dependent  Schrodinger  equation,  and  to  reconstruct  the 
pos t -CO  1 1 i s i on  wave  function  by  adding  the  scattered  packets.  If 
the  width  of  each  packet  is  narrower  than  the  length  scale  over 
which  the  potential  varies,  the  center  of  the  packet  moves 
according  to  classical  mechanics.  The  calculation  of  the 
scattering  of  one  packet  is  thus  reduced  to  the  calculation  of  a 
classical  trajectory,  plus  four  other  first  order  differential 
equations  giving  the  complex  width  and  the  phase  of  the  packet. 

In  what  follows  we  call  this  the  simplest  Heller  method,  or  SH.M. 
Even  if  hundreds  of  packets  are  needed  to  construct  the  wave 
function  SHM  requires  less  computer  time  than  the  traditional 
quantum  methods. 

While  the  resulting  theory  has  a  strong  resemblance  with 
classical  mechanics  it  contains  all  the  important  quantum 
effects;  interference  -  which  is  crucial  in  a  study  of 
diffraction  -  is  preserved  since  the  pos t -co  1 1  i  s  i  on  packet 
amplitudes  are  added  and  then  squared;  Heisenberg  principle  is 
intact  since  the  packets  have  a  width  so  that  both  the  position 
and  the  momentum  are  allowed  to  fluctuate:  finally  the 
observables  are  computed  by’  the  rules  of  quantum  mechanics  (i.e. 
by  using  operators  and  matrix  elements,  rather  than  the  classical 
definitions).  Thus  the  theory  must  be  viewed  as  a  fully  quantum 
theory  designed  to  resemble  classical  mechanics  and  acquire  thus 
the  computational  advantages  of  the  latter. 

Such  a  theory  holds  great  promise  for  several  reasons.  (a) 
There  is  no  conceptual  difficulty  in  extending  the  theory  to 
partly  disordered  surfaces^  and  dealing  with  scattering  by 


surfaces  with  short  correlation  length  will  not  strain  the 
capability  of  existing  computers.  By  comparison  a  traditional 

3 

quantum  scattering  method  encounters  severe  difficulties  when 
translational  symmetry  is  lost.  (b)  The  theory  is  easily 
extended  to  the  case  when  some  degrees  of  freedom  are  to  be 
treated  classically,  while  the  translation  of  the  projectile  is 

4 

quantized.  This  is  very  important  since  it  permits  us  to 
calculate,  n on - pe r t u r ba t i ve 1 y .  how  the  thermal  motion  of  the 
lattice  affects  diffraction.^  Once  such  developments  are 
completed  and  tested,  they  will  augment  substantially  the 
usefulnessof  atom  diffraction  as  a  sensitive  surface  probe. 

While  the  Dr o 1 s hagen-He 1 1 e r  papers  justify  high  hopes, 
their  use  of  the  simplest  version  of  Heller's  theory  is  cause  for 
some  concern.  As  is  well  known  the  signal  produced  by 
interference  effects  depends  on  the  degree  of  coherence  of  both 
the  probing  beam  and  the  target:  scattering  of  a  partially 
coherent  beam  from  a  perfectly  ordered  surface  gives  broad  peaxs. 
as  if  the  surface  is  disordered.  Intuitively  it  is  clear  that 
the  coherence  in  the  direction  parallel  to  the  surface  is  more 
important  than  that  in  the  direction  perpendicular  to  it.  since 
that  is  the  plane  in  which  the  target  is  ordered  The 
Drolshagen-Heller  theory  does  a  good  job  in  this  regard,  since 
they  use  a  sura  of  wave  packets  to  build  a  piece  of  a  two 
dimensional  planar  wave  (i.e.  full  coherence)  over  the  unit 
lattice,  scatter  it  from  the  unit  cell  and  then  construct  the 
final  wave  function  from  it  by  using.  Bloch's  theorem.  However, 
the  beam  is  partially  coherent  in  the  z-direction  since  it 
resembles  a  pl*anar  wave  only  on  a  length  scale  equal  to  the  width 
of  the  packet  (which  plays  the  role  of  a  coherence  length) 
Existing  experience  with  interference  suggests  that  if  this  width 
does  not  exceed  some  of  the  important  length  scales  of  the 
target,  the  partial  coherence  of  the  beam  will  broaden  the 
diffraction  peak;  the  use  of  wave  packets,  would  then  affect  the 
diffraction  results.  One  cannot  avoid  this  difficulty  by  using 
broad  packets  since  the  dynamic  scheme  used  for  packet 


propagation  (i.e.  SHM)  requires  them  to  be  narrower  than  the 
length  over  which  the  potential  changes. 

The  purpose  of  this  article  is  to  examine  more  closely  the 
conflicting  demands  imposed  by  the  need  for  adequate  coherence 
(i.e.  broad  packets)  and  for  computational  simplicity  (i.e. 
narrow  packets).  We  do  this  by  comparing  the  analytic 
expressions  obtained  for  the  diffraction  of  partially  coherent 
(wave  packets)  and  fully  coherent  (plane  waves)  neutron  beams  by 
a  bulk  lattice,  a  lattice  of  ordered  scattering  centers  located 
on  a  hard  wall,  and  a  set  of  scattering  centers  located  on  a  wall 
and  affected  by  static  or  thermal  disorder.  In  all  cases  we  find 
that  the  use  of  wave  packets  can  be  justified  if  their  width 
exceeds  certain  length  scales  describing  the  properties  of  the 
target,  such  as  lattice  spacing,  potential  range  and  coherence 
length  (i.e.  the  decay  length  of  the  correlations  in  the  target 
density  fluctuations). 

However,  since  neutron  scattering  can  only  provide,  at 
best,  qualitative  guidance  for  what  would  happen  in  the  case  of 
atom  scattering,  we  also  undertake  a  numerical  study  of  He 
diffraction  by  using  the  SH.M  (as  applied  to  diffraction  by 

Drolshagen  and  Heller).  We  find  that  no  matter  how  we  chose  the 

initial  width  of  the  packet  the  interaction  with  the  surface 
broadens  the  packet  to  the  extent  that  the  requirements  for  the 
validity  of  the  SHM  propagation  equations  are  not  satisfied.  A 
detailed  comparison  of  the  results  obtained  by  the  DH  method  with 
those  given  by  the  usual  quantum  theory  shows  a  good  agreement 

between  the  two.^  This  seems  to  suggest  that  SHM  works  well--at 

least  in  the  case  examined--outside  its  stated  range  of  validity. 
While  this  is  fortunate  one  cannot  help  being  slightly 
uncomfortable  in  regards  to  the  reliability  of  the  results  for 
situations  in  which  "exact’  quantum  calculations  are  not 
available  lor  obtainable)  for  testing  those  produced  by  SHM. 
Methods  that  increase  the  reliability  of  wave  packet  calculations 
are  available  and  their  application  to  diffraction  will  be 
iiscussed  in  a  future  article. 


o 


r  r  .  AN  ANALYSIS  OF  WAVE  .PACKET  DIFFRACTION  BV  L’SING  SIMPLE 

ANALYTICALLY  SOLVABLE  MODELS. 

II. A  Introductory  remarks. 

We  examine  several  simple  models  to  clarify  how  the 
diffraction  peaks  are  affected  by  the  use  of  particle  states 
which  are  defined  to  maintain  a  Gaussian  dependence  on  the  z 
variable  throughout  the  collision  process.  We  find  that  this 
practice  modifies  the  intensity  of  the  diffraction  peaks  in  a 
manner  which  depends  on  the  relationship  between  various  length 
scale-s:  the  width  of  the  packet,  the  lattice  spacing,  the  range 

of  the  potential  and  the  correlation  length  of  the  scatterers. 

After  describing  briefly  in  Section  II. B  the  model  for  the 
incident  wave  function  and  its  propagation,  we  discuss  in  Section 
II. C  elastic  neutron  scattering  by  a  perfectly  ordered  bulk 
crystal.  For  this  case  the  part i c  1  e- 1  a 1 1  i  c e  interaction  can  be 
adequately  described  by  a  potential  consisting  of  sums  of  6- 
functions  and  the  scattering  process  *  treated  by  first  order 
perturbation  theory..  The  only  lengths  in  the  problem  are  the 
width  of  the  packet,  the  de  Broglie  wavelength,  and  the  lattice 
constant.  In  Section  II  D  we  consider  a  disordered  lattice 
(static  or  thermal  disorderl.  which  adds  two  new  lengths  to  the 
problem;  the  correlation  length  and  the  mean  amplitude  of  the 
deviations  of  the  atom  positions  from  the  perfect  lattice.  In 
II. E  we  consider  a  model  that  simulates  surface  scattering  by 
looking  at  a  hard  wall  on  which  we  place  an  ordered,  two 
dimensional  array  of  scattering  centers.  We  allow  the  projectile 
to  interact  with  these  centers  through  an  exponentially  repulsive 
potential  to  study  how  the  relationship  between  the  coherence 
length  of  the  packet  and  the  range  of  the  potential  affects 
diffraction. 


I  I  .  B .  The  Model 


The  initial  wave  function,  at  a  time  t^  prior  to  the 
coilision.  is 

-oa 

4/ .  ( ? )  =  (4TTa^)  r  dK  exp{-(K-)c  )^4a^}exp{-icz  )exp<il<'?} 

1  J  I  11 

-00 

2  2 

=  exp{-a  (z-z  }exp(iK-R}  exp{iK^(z-z^)}  12.11 

Here  K  =  (k  .k  )  and  R  =  (x.y)  are  two  dimensional  vectors 
x  y 

parallel  to  the  surface,  k  is  the  z-component  of  the  wave  vector 
and  k  =  (K.k)  is  the  total,  three-dimensional  wave  vector. 

The  z  dependence  o  f  is  a  Gaussian  centered  at  z,  with  a 

-1  ‘  ^ 
width  cr=a  This  width  defines  the  spatial  scale  on  which  the 

variation  of  with  z  (at  fixed  x  and  y)  resembles  a  planar 
i 

wave;  thus  cr  is  the  coherence  length  of  the  beam  along  the  z 

direction.  To  simplify  matters  we  have  assumed  perfect  coherence 

:n  the  x  and  y  direction  so  that  cr  is  the  only  coherence  length 

in  the  problem.  I  f  a  0  .  ij)  becomes  a  planar  wave  and  the 

i 

coherence  length  becomes  infinite. 

The  initial  state  41  is  propagated  with  the  Hamiltonian 

i 

H  =  H  *V=-(h^.''2m)V^-V(r)  (2.2) 

0 

to  some  final  time  t^  chosen  such  that  at  t^  the  packet  has 
already  passed  through  the  lattice  and  stopped  interacting  with 
it  We  want  to  compute  the  probability 

=  Sd'r  exp{-il<'-r’}  Ult^-t^l  il).!"?)  ^  (2.3) 

tnat  the  scattered  particle  emerges  from  the  solid  with  a  wave 
vector  The  propagator  Ult^-tj^l  is  computed  by  using  time 

dependent  perturbation  theory;* 


U  (  t  ) 
0 


(2.5) 


=  exp{-iH  t/h>. 
o 

and  V  is  the  pa r t i c 1 e - su r f a ce  interaction  energy. 

Since  the  wave  packet  is  not  a  stationary  state  of  H  we 

o 

use  time  dependent  perturbation  theory  rather  than  tlie  Golden 

g 

Rule  expression  which  is  normally  used  in  the  theory  of  neutron 
scattering . 
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I  I  .  C  .  Neutron  scattering  bv  a  three-dimensional  lattice. 

Neutron  scattering  can  be  adequately  described  by  using  the 
n eu t r on  -  1  a 1 1 i c e  interaction 


n  ,  in  .  1  =  1 


where  the  vectors  r  ,  =  c(nx-my-lz)  give  the  positions  of  the 

lattice  atoms,  c  is  the  lattice  constant,  x.y.z  are  unit  vectors 
along  the  coordinate  axis  (we  consider  a  cubic  lattice  with  the 
crystal  axis  along  x.  y.  z).  n.m.l  are  integers.  N^,  N^.  N^  are 

the  number  of  atoms  in  the  x.y  and  z  directions,  and  N  is  the 
total  number  of  lattice  atoms.  The  crystal  is  located  between 
z=0  and  and  is  infinite  in  the  x  and  y  directions.  Since 

we  take  in  Eq.(2.1)  z^<0  and  k^>0  the  wave  packet  starts  in  the 
half  space  z<0  and  moves  towards  the  z=0  face  of  the  crystal  with 
the  group  velocity  v^  =  fiK^/m. 

Using  the  Eqs.  (2. 1-2.6)  we  can  calculate  the  probability 
PtiT  )  that  a  neutron  whose  initial  state  is  described  by  Eq 
(2.1)  will  have,  after  the  collision  with  the  lattice,  a  final 
momentum  k’  The  result  is 

P(l<')  =  TTIV  hav,)^  f(Ak  )  f(Ak  )  F  I2.T) 

o  1  X  y 

with 

f'Ak  )  =  (sin(cik  N  2!''N  sinlcAk  '  2  ]  .  'ri  =  x  or  v.  (2.3! 

a  iu  4  (u  la 

and 

F  =  N"“exp(-(l,2)(Aehav  )^}  I  Z  {exp[-iAic(l-l')c] 
z  ^11' 

•exp[iA£(l-l')c  nv  ^  j  ) 


12.9) 


c  -< 


1 


is  the  (iiffere.Tce  between  the  ?.  -  components  of  the 


final  momentum  and  the  initial  mean  momentum  of  the  packet,  hk 


k  -k  .  ^  =  .x.y,  is  the  momentum  transfer  in  the  x  or  y  direction; 

t  O')  2  ^  ">  5 

and  Ac  =  (n  2m)(k'  -kl.  with  k’  =  (k  )  -  (k  i”  -  <  '  and  k"  = 

0  0  0  1  X  y  1 


k'-k'-<",  is  the  ener^iy  transfer.  Since  we  do  perturbation 


theory  the  packet  maintains  the  velocity  v,  as  it  'Joes  throug'n 
the  lattice.  In  carrying  out  the  integrals  over  time  we  have 

assumed  that  the  width  of  the  packet  remains  constant  in  time. 
This  is  not  quite  correct  since,  in  perturbation  theory,  tne 


packet  width  has  its  free  space  variation  given  by  a(tl  = 


,  2  -  1  '  ’ 

1-2 ha  it  m]  “.  The  neglect  of  this  time  dependence  allows  us 


to  get  a  simple  result,  hopefully  without  altering  substantia 
the  physics  of  the  problem. 


In  order  to  better  understand  the  modifications  introduced 
by  the  use  of  a  wave  packet  we  compare  Eqs  (2.7-2.91  to  tli e 
diffraction  formula  for  neutrons  whose  initial  state  is  a  planar 
wave.  The  transition  rate  for  the  latter  case  is 


W(k'  ) 

= 

V“  {  2-T 

0 

h")  f(Ak  ) 

X 

f(Ak  )  f 

y 

(  Ak 

2 

)  6  (  u  (  k  ' 

)  -Ui  ( 

k  !  ! 

i  2  .  10) 

s  N  ,  N  ,  N 
\  y 

2 

O 

O 

infinity  f  (  Ak  '■  -*5  (  k 

-  k 

-  P 

-  2Trn  c  ! 

w  t;  e  r 

e 

n  =  0  .  1  .  =  2  , 

and 

u  =  X . y  or 

2 ;  thus  . 

the 

p  r  o  d  u  c  t 

f  (  A  k  .  i  r  (  A  k 
s 

y 

)  f  (  A  k  ) 
2 

genera  tes 

Bragg's 

di  f  f 

'.'a c  '  ion 

c  0  n  d 

i  t  :  on  . 

F'.irthermoro 

since 

we  conside 

r  e  I  a  s  t  i 

c  scattering 

.  d 

f  -  f  u  n  c  t  '.on 

enforcing 

energy  conservation 

appears 

i  n 

Eq . ( 2 . 1 0 

I 

The  essential  difference  between  a  wave  packet  and  u  pl.anar 
wave  initial  state  is  that  the  term  h  =  6(c'-c^-2TTn  cloiAci 
appearing  in  the  planar  state  case  and  representing  t(\e 
conservation  of  the  2-component  of  momentum  and  of  the  energy  is 
r-^placed  by  F  Since  the  wave  packet  is  not  an  eigenvalue  of  the 
.tern  order  Hamiltonian,  or  of  the  .momentum  operator,  it  is  nor 
3 u  r  p  r  .  s  1  n g  ‘hat  the  simple  conservation  conditions  e  .x p  r  e  s  s e d  by 


h.  for  planar  wave  initial  condition,  are  no  longer  valid.  .Vote 
that  as  a  -*  oo,  F  -*•  h  .  as  it  should. 

The  manner  in  which  the  use  of  a  wave  packet  affects 
diffraction  can  be  understood  by  e-Kamining  F.  Carrying  out  the 
double  sum  and  taking  -►  «  we  obtain 

F  =  e.xp  { -2  (  Ae /fiav  )  ^  )  6  (  Ak  -  2TTn/ c -Ae  hv  ^  )  (2.11) 

In  the  case  of  a  planar  wave  the  permissible  final  momenta  can  be 
found  by  specifying  n^  and  n^  since  the  restriction  imposed  by 
the  energy  conservation  term  6(Ae)  fixes  n^  and  thus  the 
permisible  values  of  k'.  Since  the  latter  restriction  disappears 
when  a  wave  packet  is  used,  a  larger  number  of  values  for  k'  are 
allowed  for  a  given  pair  {n^.n  ).  Thus,  what  used  to  be  a  sharp 
Bragg  spot,  when  the  initial  state  was  a  planar  wave,  may  degrade 
into  a  large  number  of  spots  having  the  same  azimuthal  angle,  but 
different  polar  angles.  If  the  detector  does  not  resolve  these 
spots  they  form  a  strip  whose  intensity  1(e)  varies  with  the 
polar  angle  e.  The  spots  are  determined  from  Eq  .  (2.11)  through 
the  condition 

Ak  =  2Trn/ c  -  Ae/hv,  n  =  (n.n)  (2.12) 

n  n  X  y 


where  n  and  n  are  fixed.  The  intensity  of  the  spot 
X  y 

corresponding  to  Ak  is  determined  by  the  magnitude  of  the 

n  2  2  '  2  2 

Gaussian  exp{-2(Ae  hav,)  ).  Since  Ae  =  (h  2m)(K  -K  ) 

,  «  1  n  n 


-  ( 2m ) ( K  ^ 

2  n 

-  k2, 

2 

.  the  -  only 

part  of  Ae 

n 

depending  on  Ak 

i  s 
n 

a 

n 

5  (h  2m)(2AK 

k*Ak 
n  n 

)  . 

The  variation  of  the  intensity 

with  Ak 

i  s 
n 

then  given  by 

e  X  p  {  - 

2‘“n 

/hav  ^ ) 

} .  Since  k 

is  of  order 

0  f 

the 

d  e 

Broglie  wave 

vector 

‘^d 

and  in 

diffraction 

studies  we 

use 

a  d  e 

Broglie  wave 

length 

0  f 

order 

of  the  lattice  spacing, 

we  f 

i  nd 

that 

ihe  exponential  cuts  off  the  values  of  Ak  which  are  larger  than 

12  1  '■  2  ^ 

'.<3;  --  (  2TTa,  c  1  The  wave  packet  gives  the  same  result  as 

a  plan.ar  wave  when  k  a  is  very  small.  This  means  that  a 


must 


■■  "I 


1 1 


be  larger-than  the  de  Broglie  wave  length  and  lattice  spacing. 

..Vote  that  besides  turning  a  spot  into  a  strip  along  the 

polar  angle,  the  use  of  a  wave  packet  also  changes  the  intensity 

of  the  spots  as  (n^.n^)  are  changed,  through  the  Gaussian  term 

appearing  in  F.  This  effect  is  negligible  for  those  peaks  for 

which  Ae  /hav,  is  very  small.  Since  Ae  is  less  than  the 
n  1  n 

incident  energy  we  have  an  upper  bound  for  Ae  leading  to  the 

2  2  - 1  5 

condition  h  k!"/ 2mha  v  ,<  <  1  which  gives  a  <<k.,k..  If  k.  is  of 

1  .1  1  _1  i 

order  of  k^.  if  a  <<  k^  the  effect  of  the  use  of  a  wave  packet 
on  the  spot  intensity  is  negligible 

The  manner  in  which  the  use  of  the  wave  packet  affects 
diffraction  can  be  understood  by  examining  what  happens  in  the 
double  sura  in  Eq.  (2.9).  If  the  exponential  exp(-iAe(l-l')c/hVj} 
was  absent,  the  double  sura  would  give-because  of  constructive  and 
destructive  i n t er f e r enc e - the  Bragg  condition  for  Ac.  The 
presenc*  of  that  exponential  alters  the  perfect  interference 
However,  since  the  packet  looks  like  a  planar  wave  over  the 
length  a  ^  this  alteration  does  not  lead  to  a  complete 
destruction  of  the  interference  pattern.  To  see  this  let  us 
divide  the  values  of  11-1')  into  two  sets.  One  set  contains 
those  values  of  1-1'  for  which  exp{-i(l-l'IAec  hv^)-!.  i.e.  those 
satisfying  (  1  -  1  '  )  A  e  c hv  <  1.  The  other  set  contains  the 
remaining  values  of  1-1'.  Thus  if  we  denote  (hv^'cAel  =  Al.  tne 
first  set  is  defined  by  l-i'5  Al.  .Vow  we  can  write  the  double 
sum  appearing  in  F  as 

Al 

Z  exp{-iAK{l-l')c}  -  remaining  terms. 

1-1 ' -0 

In  the  remaining  terms  the  phase  Ae(l-l')chv^  is  larger  than  2  tt 
and  thus  acts  as  a  random  phase  leading  to  cancellations  in  the 
double  sum  and  destroying  interference.  Thus  interference  is 
present  in  the  first  term  only  and  we  have 


.V 

z 

ot  Z 
1=0 


2 

a  S  Z  exp{-iAK(l-l')c}exp{-2(^iehav  )  } 

1=0  1-1  ' =0 

-  background  . 

2 

The  double  sum  is  proportional  to  {sin(AKcAl/2),Alsin(AKc2)} 
which  has  a  central  peak  at  Ak=0.  having  a  width  A1  ^  .  Now  since 

Ai-hv  /cAe  and  the  Gaussian  cuts  Ae  off  at  hav  ,  the  width  of  the 

^  - 1  ^ 
function  FIAk)  is  (Al)  --  chav^/hVj^  =  ca.  If  ca  is  very  small 

li  e.  if  O'  c)  the  packet  scattering  resembles  that  of  a  planar 

wave  . 


This  example  was  used  here  only  because  it  is  the  simplest 
possible  situation  in  which  we  can  examine  the-  effect  of  using 
wave  packets  on  diffraction.  Since  in  atom  surface  scattering, 
which  is  the  phenomenon  in  which  we  are  interested,  the  atom  does 
not  pass  through  the  solid,  the  effects  discussed  above  are  not 
directly  relevant  to  atom  surface  scattering.  However,  the 
mathematical  machinery  through  which  these  effects  are  generated 
is  the  same.  A  modei  which  mimics  the  physical  conditions  of 


surface  scattering  is  discussed  in  Section  II. E. 
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It  is  intuitively  clear  that  the  presence  of  partial 

disorder  in  the  position  of  the  scattering  centers  must  affect 

the  diffraction  pattern  in  a  manner  similar  to  the  partial 

coherence  of  the  beam.  Since  one  of  the  hopes  opened  up  by  the 

Drolshagen-Heller  paper  is  the  possibility  of  interpreting  atom 

1 

scattering  from  partly  ordered  surfaces',  it  is  important  to  use 
simple  models  in  order  to  try  to  understand  how  the  use  of  wave 
packets  might  influence  the  results  of  such  a  computation.  The 
key  practical  matter  is  to  find  under  what  conditions  the  use  of 
a  wave  packet  for  the  2  dependence  of  the  initial  state  leads  to 
the  same  results  as  the  use  of  a  planar  wave. 


We  consider  here  lattice  disorder  that  can  be  described  by 
a  small  displacement  I'  (in  what  follows  we  use  Bsin.m.ll  and 

O 

8'=(n’,m'.l'))  of  each  solid  atom  relative  to  the  lattice  point 
.  (i.e.  the  position  of  the  atom  located  near  the  site  B  is 

U  O 

=  t-U).  Such  disorder  could  be  caused  either  by  thermal  motion 
o  5 

or  by  poor  crystallization  when  the  solid  was  formed  li.e.  static 
disorder).  We  assume  that  in  both  cases  is  a  Gaussian 
variable,  and  therefore  it  can  be  completely  characterized  by  the 
correlation  function  D(B-B';  tl  =  <L'  (t)  L’  .(01>  and  the  condition 

D  O 

<l'(t;>  =  0.  In  the  case  of  thermal  disorder  U„  is  a  Gaussian 

O  o 

variable  if  the  solid  is  harmonic.  For  static  disorder  (:h  which 
case  the  displacement  is  time  independent)  we  pos.culate  the 
Gaussian  property;  the  true  statistical  properties  of  depend 
on  the  mechanism  of  crystallization. 


The  disorder  introduces  two  new  length  scales:  the 

1  '  ■’ 

amplitude  A=<U  'J.>  '■  and  the  correlation  length  1  defined  as 

D  D  C 

the  distance  beyond  which  the  density -  density  correlation 
function  e'.^uals  the  density  squared  ii  e.  <.oir)DlO)>-<pir!><otO)> 
for  r  >1^).  Here  <>  is  either  a  thermal  average  or  an  average 


over  the  distribution  function  characterizing  the  static 
disorder . 


To  understand  how  the  use  of  a  wave  packet  in  the  z 
variable  affects  the  diffraction  by  a  disordered  lattice  we 
consider  neutron  scattering  by  a  bulk  crystal.  We  use  the 
potential  defined  by  Eq  .  ( 2 . 6  )  and  write 


(2.18) 


Using  the  methods  described  earlier  in  this  article  combined  with 

g 

standard  neutron  scattering  methods  we  find 


PtiT'.ui)  =  (V  /hN)^  J*  dt  J"  dt'  exp{iAe(t-t')  h). 

0  o 

(2.19a) 

Z  Z  exp  {- iAl<«  (T- S  (  B  ,  B  ;t,t') 

B  B'  ®  ® 

with 

S(B.B';t-t')  =  <exp{ifil<-(J_,(t’)}  exp{-a^(l'c-Un.(t’)*z- 

o  o 

z(t'))^}  exp { - iAk • U- ( t  )  }  exp{-a^(lC'U„(t)'Z-z(t)  )~> 

O  D 


(  2  .  19b  1 


Here  iik=k'-k,.  z  (  t )  =  z  *hK  ,  t  /  m  and  U_(t)-z  is  the  projection  of 

tiie  displacement  vector  at  site  B  along  the  z  direction.  The 

type  of  average  depends  on  whether  we  deal  with  static  or  thermal 

disorder.  If  we  use  a  planar  wave  in  the  z  direction  S  becomes 

the  structure  factor  of  the  lattice  (essentially  the  density- 

8  9 

density  correlation  factor).  '  Note  that  because  of  the  use  of 

the  Gaussian  wave  packet,  the  displacement  amplitude  U  appeares 

o 

squared  in  some  of  the  exponents  present  between  the  averaging 
brackets.  This  makes  the  evaluation  of  the  averages  cumbersome. 
For  this  reason  we  neglect  the  terms  containing  the  square  of  the 
displacement,  and  then  compute  the  average  by  standard  methods. 

It  can  be  shown  that  this  does  not  introduce  substantial  errors 


I 


We  can  thus  compute  S(B.B  ;  t.t’)  and  obtain 

S(B.B';t.t')  =  exp{ -a^[  (  Ic-z  (  t )  )  ^- (  1  '  c-2  (  t  ’  )  )  }  ■ 

exp  {  [  -  2  {  Ak  )^-4a^(l  ’c-zCt'  )  )^-4a^(  lc-z(  t)  )“-4il  )  a^ 

(  (  1-1  ’  )c-z(  t  ’  )-z( t)  )  ]  D(B.B’:0)/2}  • 

exp{  [2(Al<)^-aa“(l  'c-z(t'  ))(  lc-z(  t)  )  - 
4i(Ak-z)a^(c(l'-l)-z(t)-z(f))]  D(B.B’:t-t':)/2l  . 

(2.20) 

In  the  case  of  thermal  disorder,  the  correlation  function  D  is 
given  by 

0(B.B-  : t)  =  <Ug(t)Ug,  (0)> 

=  (fi/12NM)  1  u('q)  ^  {  e  .x  p  [  i  q  ■  ( T  ^  -  ?g)-iu)(q)t) 

q 

(  n  (  q  )  -  1  ) -exp  C  -  iq  •  (Tg  , -r’g  )  -  iul  q  )  1 1  n  (  q  )  >  .  (2.21) 

where  q  and  u(q)  are  the  phonon  wave  vector  and  frequency, 
respectively.  .M  is  the  mass  of  the  lattice  atom  and  n(q)  = 
[l-exp(-Fiu(q),kT)]  ^  is  the  thermal  phonon  population. 

The  exponential  in  £q.(2.20)  containing  D(B.3';0)  is  the 

analog  of  the  Debye-Waller  (DW)  factor.  If  a  planar  wave  is 

used,  instead  of  a  wave  packet,  that  term  reduces  to  the  usual  DW 

-*■  2 

factor  exp(-2W)  =  exp{-2(Ak)  D(B.B;0)}  (we  assume  a  simple 

Bravais  lattice).  The  exponential,  in  Eq .  (2.20).  containing 

D(3.B':t-t')  represents  the  effect  of  the  inelastic  collisions. 

For  Che  purpose  of  analyzing  what  happens  to  the  DW  factor  when  a 

wave  packet  is  used  we  can  consider  the  elastic  term  only:  this 
3  9 

:s  done  by  taking  D(B.B';t)=0. 


We  can  now  compu'te  P(k')  for  an  elastic  collision  and 

obtain 

Pik'  )  =  TTIV  hav  ,  )  ~  (  1 -4a"D  )  '  ‘  ex  p  {  -  2  W  (  1  -  4  a  '  D  )  '  *  ! 

0  1 

exp(-[2(l-4a^D)l  ^(Aenav^l^} 

exp{2  Ak  (  ie  )  D  [  V  .ln(  1 -4a^D  )  ]  '  *  )  fiAk  )f(Ak  )  (2.22! 

-  X  y 

N  "  I  exp{-i[K'-<)-(Ae  hv,)]{l-l)  c 
z  1  1  ,  1 

The  above  equation  is  valid,  due  to  some  appro xi. nations  we  made. 

.j  n  1  y  if 

4a"D(B.B;0)<l.  (2.23) 

We  see  that  the  use  of  a  packet  results  in  the  replacement  of  the 
customary  DW  factor  exp(-2W)  with  exp[-2W(l-4a  D)  ].  Thus  the 
partial  coherence  in  the  beam  leads  to  greater  reduction  of  the 
intensity.  Intuitively  this  is  not  surprising  since  we  can  thinx 
of  the  scattering  of  a  planar  wave  by  a  perfectly  ordered  solid, 
as  the  scattering  of  a  fully  coherent  probe  by  a  fully  coherent 
object.  Introducing  disorder  in  the  object  is  equivalent  to 
turning  it  into  a  partially  coherent  target.  The  use  of  a  w.o'.-“ 
packet  causes  some  partial  destruction  of  the  probe  ;ohe:-ence. 
which  is  equivalent  to  a  "disorder"  in  the  beam.  We  expect  rha: 
the  scattering  event  will  convolute  the  incoherence  of  the  b'-am 
with  that  of  the  target;  thus  the  intensity  decrease  associated 
with  the  lattice  disorder  is  accentuated  by  the  ’disorder'  of  the 
beam.  If  4a"D(B.B:0)  approaches  1.  the  scattering  probability  is 
strongly  diminished.  This  indicates  that  unless  A  the 

Debye-Waller  factor  obtained  by  a  calculation  using  a  wave  packet 
:n  the  z-direction  is  substantially  larger  than  the  one  obtained 
■.'H'-i  a  planar  wave  is  used.  If  the  lattice  temperature  is 
several  hundred  degrees'Kelvin  a  beam  coherence  length  of  only 
in-  A  will  lead  to  noticeable  changes  in  the  Debye -Waller  fa^t'r 
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Since  the  form  of  the  potential  is  such  that  SHM  propagation 
requires  a  packet  width  of  about  1/4  A.  the  use  of  the  wave 
packets  could  lead  to  either  errors  in  propagation  (when  a  large 
width  is  used)  or  in  the  DW  factor  (when  a  small  width  is  used). 
We  will  come  back  to  this  in  the  next  section. 

As  in  the  zero  temperature  case  we  get  in  Eq  .  (2.22)  a 

Gaussian  in  Ae  reflecting  the  effect  of  the  energy  spread  in  the 
initial  wave  function  the  width  of  this  Gaussian  has  been 

however  modified  by  disorder  and  decreases  as  the  fluctuation 
amplitude  A  increases.  There  is  also  a  new  term  which  raises  or 
lowers  the  intens'ty  of  various  peaks  depending  upon  the  sign  of 
Ak  and  Ae .  For  small  a.  Ae  is  almost  0  and  this  term  is  small 
compared  to  the  one  modifying  the  Debye- Waller  factor.  The  Bragg 
conditions  are  modified  as  in  the  zero  temperature  case. 

In  summary,  the  use  of  a  wave  packet  in  the  z-direction  has 
the  effect  of  increasing  the  Debye  Waller  factor,  and  generating 
a  new  term  which  further  modifies  the  amplitudes  for  Ae?^0.  As  in 
the  zero  temperature  case  the  restrictions  concerning 
conservation  of  energy  and  z-momentum  are  loosened,  resulting  in 
peak  broadening  along  the  polar  angle  ©.  The  peaks  are  centered 
at  positions  specified  by  e(K’)  =  sic). 

Before  proceeding  further  we  make  a  few  remarks  concerning 
some  of  the  practical  problems  encountered  in  computations.  We 
have  so  far  concentrated  on  the  role  of  the  partial  coherence 
along  the  z-direction;  however,  when  numerical  calculations  are 
carried  out.  some  difficulties  appear  in  treating  fully  the 
parallel  (i.e..  x.y  direction)  coherence  if  the  lattice  is 
disordered  or  it  undergoes  thermal  motion.  For  a  perfectly 
ordered  target  we  can  build  up  the  wave  function  of  the  scattered 
neutrons  by  writing  the  "parallel"  planar  wave  exp(iK-R)  as  a  sura 
of  wave  packets  and  propagating  only  those  packets  that  hit  a 
given  unit  cell:  Bloch's  theorem  allows  us  to  construct  the 


n  f  ?  ft  ? 


T  V>  o  c  m  o 


nrororinro  ran 


‘  ■  I  '.  :  -^  .  1  .  I  i  a  '  :  n  •> :  'he  .;  e  n  y  r  -  ^  .  r  .•  r  ‘  ■ 

A’  p  r  r  r  .,1 1' ba  t:  i  on  theory)  the  soiiai-.e  of  the  “  o  u  r  .  e  ;- 

of  cue  expectation  value  of  Che  density  operat or 


"♦  -qH-1 

<o(r)>  =  Trie  Z  Z  olr-r.)} 

B  ® 


One  can  easily  show  that  <p(r)>  is  a  periodic  function  consistin 
of  sums  of  Gaussains  centered  around  Che  lattice  points. 
Therefore  one  can  construct  the  Debye- Waller  effect  by 
calculating  scattering  by  a  unit  cell  whose  atoms  are  allowed 
move  thermally,  and  using  Bloch's  theorem.  If  we  select  only 
elastic  events  this  calculation  will  give  the  effect  of  the 
temperature  on  Che  peak  intensity,  but  not  on  its  width. 


Some  difficulty  appears  when  inelastic  effects  are  inciu 

to  compute  the  effect  of  the  temperature  on  the  width  of  the 

diffraction  peaks.  In  a  simple  perturbation  theory  such  effec 

appear  in  the  cross  section  through  the  Fourier  transform  of 

'.  pl?,t)p(0.0)>  with  respect  to  time  and  space.  The  projectiie 

samples  now  the  correlated  motion  of  all  t  'he  p  a  r  t  i  c'  1  e  s  w  i  *  h  :  n 

density -density  correlation  length  1  .  Therefore  we  must  sc. at 

^  •) 

a  two  dimensional  planar  wave  from  a  surface  of  area  1'  and  th 

retjuires  the  propagation  of  a  very  large  number  of  Gaussians. 

simultaneously  the  generation  of  the  correlated  motion  of  all 

0 

1  1 1  i  c  e  .atoms  within  1  '  .  This  is  -a  rather  e  x  u  e  n  s  i  v  e  c  a  1  c  a  i  a  t  : 

c 

which  will  provide  a  great  incentive  for  finding  a  way  i  i 
circumvent  it  by  incorporating  the  inelastic  processes  throug.c 
P'.'rturbation  theory:  this  should  work  for  light  scatterers  and 
ri  e  a  V  y  !  a  t  t  ; 


c  e  atoms. 


[  I  .  E  .  Surface  scattering 


In  this  section  we  examine  the  effect  of  the  beam  coherence 
<j  on  scattering  from  a  surface  consisting  of  a  hard  wall  on  which 
we  superimpose  a  potential  V(x.y,z)  which  is  periodic  in  a 
direction  parallel  to  the  surface  and  decays  exponentially  in  the 
direction  of  the  vacuum.  To  compute  the  scattering  probability 
we  use  the  distorted  wave  approximation;  the  hard  wall  is 
treated  exactly  and  V(x.y,z)  is  considered  a  small  perturbation. 
This  model  allows  us  to  gain  some  insight  into  the  relationship 
which  the  coherence  length  in  the  z  direction  ct  must  have  with 
the  potential  range  X  ^  and  the  lattice  constant  c.  in  order  that 
scattering  of  the  wave  packet  by  the  surface  mimics  the  case  of  a 
planar  wave.  We  find  that  unless  a  exceeds  both  X  ^  and  c  the 
calculated  diffraction  intensity  is  strongly  affected  by  the  use 
of  a  wave  packet  . 


A  zeroth  order  Hamiltonian  which  describes  surface 
scattering  is  obtained  by  adding  a  term  V'(z)  to  in  Equation 
(2.2).  where 

V”  (  z  )  =  «  z  <  0 

and  (2.12) 

V '  (  z  )  =  0  z  >  0  . 

The  eigenstates  of  H  *  V'(z)  are 


41^  ("r  I 


i KR  ,  iKZ 
e  (  e 


-  iKZ  , 

e  ) 


z  >  0 


and 

4t  .  (  T )  =  0 
1 


z  <  0 


12.13) 


The  initial  state  is 
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tlJ.(r)  =  (4TTa^)  X  dK  exp{-[<ic-K^)/2a]"} 


-i(K  -zl  iK-R  iKz  -iKz.  „ 

e  1  le  (e  -e  ),  forz>0 


(2  14) 


and 


»|i  .  (  r  )  =  0  . 


for  z  i  0 


The  wave  function  if  ^  (  r )  consists  of  two  wave  packets,  one 
in  the  vacuum  and  the  other  one  inside  the  solid,  symmetrically 
located  with  respect  to  the  surface.  When  lii  .  (  r  )  is  propagated  in 
time,  both  packets  move  towards  the  surface  and  overlap  there. 
Then  the  packets  separate  again  and  there  is  only  one  in  the 
vacuum.  Since  by  definition.  tl>^(r^)=0  for  z  S  0.  the  sum  of  the 
two  packets  has  no  meaning  for  z  S  0;  only  the  values  of  41.  for  z 
>  0  contribute  to  the  wave  function  of  the  system. 


To  the  zeroth  orde'r  Hamiltonian  described  above  we  add 


V(  r  ) 


n  ,  m 


V 

xy 


(  r-R  ) 
nra 


where  R  =  nc.x  *  mcy 
nm 
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This  potential  introduces  an  interaction  length  y  ^  in  the  z- 
direction  and  a  lattice  parameter  c  in  the  .x.y  plane. 


By  using  the  methods  described  in  the  preceding  sections  we 

find 


P  (  k  ■  ) 


4.2  2  2 
h  a  V  ^ 


V  (AK):f(Ak  )  f(Ak  )exp{-(Ae/ahv)2} 
xy  X  y  1 

(2.16) 


4{k^  -  e) 


where  e=Ae  liv,,  and  V  (AK)  is  the  Fourier  transform  of  the  x.y 
1  xy 

part  of  the  potential.  Note  that  the  Laue  conditions  for 

scattering  in  k  and  k  ,  i.e.  f(Ak  )f(Ak  ).  factor  out.  This 
X  y  X  y 

happens  because  we  have  separated  the  x,y  part  of  the  potential 

from  the  z  part  (i.e.  because  v(x.y.z)  has  the  form  f(x,y)  g(z)l 
and  because  iji.ir)  is  a  plane  wave  in  the  direction  parallel  to 
the  surface.  However,  similarly  to  the  case  of  bulk  scattering, 
the  use  of  a  wave  packet  in  the  variable  z  causes  a  broadening  o 
the  diffraction  peaks  in  the  variable  k'  (which  is  equivalent  to 
Che  polar  angle  e). 


For  the  case  of  an  infinite  surface  we  can  rewrite  Eq 
(2.16)  as  a  sum  over  diffraction  peaks 
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In  the  coherent  limit,  a  goes  to  0  .nd  the  Gaussian  term  in  Eq 


22 


(217)  dominates  the  diffraction  structure  alon?  <:  .  becoming  a 

delta  function  which  enforces  energy  conservation  for  each  peak 

The  probability  peaks  sharply  at  e  =0,  and  the  Lorentzian 

n  m 

term,  which  results  from  the  Fourier  transform  of  the  e.-tponential 
repulsion,  partially  determines  the  intensity  at  each  peak. 


For  the  case  of  a  finite  coherence  length,  both  the 
Gaussian  and  the  Lorentzian  terms  in  (2.17)  can  contribute  to  the 
diffraction  structure.  Which  term  dominates,  and  under  what 
conditions  we  reach  the  coherent  limit  depends  not  only  upon  the 
packet  width  a  ^  and  the  interaction  length  y  but  also  on  A 

n  m 

and  The  requirements  for  reaching  this  limit  (where  we  have 

a  single  peak  at  e  =0  for  each  n.m)  are  in  general 

nm 

c  o  mp  1  i  c  a  t  e'd  .  but  a  few  observations  can  be  made. 


The  Lorentzian  terra,  which  has  a  width  y ,  peaks  when  the 
perpendicular  final  momentum  c  satisfies. 


< 


nm 


We  have  assumed  that  the  uncertainty  in  k  .  which  is  equal  to 

0x2  =  is  less  than  k,.  such  that  (c  <k  .  and)  only  the 

1  nm  I 

second  term  in  the  denominator  of  Eq.  (2,17)  gives  peaks.  The 
Gaussian  term  has  a  width  >12  a  =  >12  'cr  and  peaks  at 


where  e  =  0.  We  see  that  these  three  peaks  only  overlao  for 
n  m . 

the  specular  case,  where  A  =  0  and  k  -  -c,.  If  the  packet 

n  ra  1 

width  is  too  narrow,  (making  a  too  large)  such  that  the  Gaussian 

term  is  broad  enough  to  overlap  significantly  with  the  Lorentzian 

peaks,  than  these  peaks  can  dominate.  In  this  case  we  would  see 

two  ma.'tima  as  we  scan  along  e.  one  on  each  side  of  the  e.xpected 

;for  planar  wave)  £  =  0  position.  Consider  the  case  of  some 

n  m 


low  order  diffraction  peak  where  A  is  small  and  positive,  such 
1/2  " 

that  A  '  <<  !k  i.  I f  7  <  !K we  recover  the  coherent  limit 

nm  1  1 

when  the  width  of  the  Gaussian  is  less  than  the  spacing  between 

1  '  2 

the  Gaussian  and  Lorentzian  peaks  .  i  .  e  .  ct  >  [2/A  1  Thus. 

)  1  '  2 

when  the  packet  width  is  larger  than  (2/A  '  we  get  a  single 

n  m 

narrow  peak  at  the  e  =  0  position.  For  larger  A  's  where  y  < 

nm  _  ^  n m 

A^^/'2  the  interaction  length  f  becomes  an  important 

factor.  In  this  limit  the  Lorentzian  peaks  are  proportional  to 
-  2 

■y  Thus,  for  a  longer  interaction  length,  we  generally  need  a 

broader  cr  so  that  the  Gaussian  is  narrow  enough  to  dominate  the 

scattering  and  produce  a  single  coherent  peak. 

For  various  parameter  values,  the  exact  requirements  vary. 

However,  for  the  typical  case  of  a  thermal  molecular  beam  hitting 

a  surface  with  c=2.5A  at  a  30°  angle,  with  y  ^ = 2A .  we  recover  the 

coherent  case  for  most  diffraction  peaks  when  ff  is  larger  than  a 

few  A  or  so.  That  is.  in  this  limit  we  only  get  a  single  peak 

along  e  center ed  at  the  k'  value  given  by  e  =  0.  For  a  larger 

n  m 

y.  or  shorter  Interaction  length,  this  limit  is  reached  for 

smaller  values  of  the  packet  width. 

N’ote.  however,  that  the  peak  heights,  at  their  e  =  0 

-2 

locations  are  proportional  to  a  and  all  other  dependence  on  a 

drops  out.  Thus  the  relative  peak  amplitudes  for  =  0  are 

independent  of  the  packet  width.  This  is  a  result  of  the 

assumption  that  a(t)  is  a  constant.  If  we  make  the  simple  but 

reasonable  assumption  that  the  incoming  and  scattered  packet  have 

different  widths,  a  ^  and  a'  ^  respectively,  then  for  e  =  0 

n  m 

P(k')o£{[y(a’-a)/a'a]^*[(K'-K^/a-(K'*iCj^)'a]^) 

(2.17) 

and  the  relative  probabilities  between  different  e  =0  peaks  is 

n  m 

a  function  of  the  packet  width.  For  a  realistic  a(t),  this 
dependence  is  in  general  very  complex.  When  a  and  a'  are 
sufficiently  large  such  that  the  Gaussia.-.  controls  the 
scattering,  this  dependence  goes  away. 


Numerical  studies  using  SHM  as  implemented  bv  Drolshage'-i 
and  Heller. 


The  analysis  of  the  neutron  scattering  -nocei  shows  that  the 
coherence  length  in  the  z-direction  must  exceed  the  lattice 
spacing,  the  de  Broglie  wave  length,  and  the  ra.nge  of  the 
potential.  Otherwise  the  predicted  diffraction  intensities 
depend  on  detailed  properties  of  the  packet  which  are  not 
physical  properties  of  the  beam. 

Such  requirements  cannot  be  satisfied  (within  SH.Ml  by  using 
e.xtremeiy  broad  wave  packets  since  this  approach  would  clash  with 
the  de.mand  made  by  the  method  used  to  propagate  the  packets  To 
describe  Che  situation  very  briefly  let  us  assume  that  ciie 
initial  wave  function  has  been  written  as  a  sum  of  Gaussians. 

SH.M  propagates  each  Gaussian  in  the  sum  independently;  to  do  this 
it  uses  the  time  dependent  Schrodinger  equation  with  the 
potential  expanded,  to  second  order,  in  a  Taylor  series  around 
the  instantaneous  position  r^  of  the  packet.  The  expansion  is 
iegicimace  only  if  the  potential  is  almost  constant  around  r.  on 
ciie  length  scale  set  by  the  packet  width  (in  the  coordinate 
representation)  .nc  the  time  t.  Thus  a  correc"  treatment  of  beam 
coherence  in  the  z  direction  requires  broad  packets  in  that 
direction,  while  the  simple  propagation  scheme  requires  them  to 
be  narrow.  The  solution  to  these  conflicting  demands  mignt  be  to 
construct  the  broad  packets  required  to  describe  the  incident 
wave  function  as  sums  of  very  narrow  wave  packets.  The  trouble 
with  this  idea  is  that  once  these  packets  interact  with  the 
surface  we  have  no  control  over  their  width.  In  particular',  a 
narrow  p.ncket  .may  be  broadened  by  the  interaction  with  ti\o 
surf. ice-more  then  the  propagation  scheme  would  safely  permit. 

The  si. m pie  .models  used  in  the  previous  section  cun  not  address 
such  issues  since  they  use  perturbation  theory  and  the  packet 
■vrdth  changes  in  time  according  to  the  free  space  dynumics. 
Therefore  to  find  out  whether  our  fears  are  founded  we  must  turn 


to  numerical  calculations. 

The  calculations  used  here  were  carried  out  like  in  the 
Orolshagen -Heller  papers  (DH)^  which  use  SHM. 

In  representing  the  planar  wave  as  a  sum  of  Gaussians  we 
use  packets  of  the  form 

glT.t)  =  e.xp{-  [  ( "r-r  )  •  A(  t  )  •  (’r-'r  ^  ) -k^  •  (T-T  ) -r  ]  } 
fl  z  t  t  t  t 

where  r^.  k^  are  real  functions  of  time  representing  the  mean 

position  and  momentum  of  the  packet:  A(t)  is  a  three  dimensional 

matrix  with  time  dependent  complex  elements:  is  a  complex 

function  of  time.  At  a  time  prior  to  the  collision  with  the 

surface  such  packets  are  added  up  to  form  a  planar  wave  function 

in  the  x  and  y  variables:  a  single  packet  is  used  in  the  z- 

variable.  Since  prior  to  the  collision  the  x.y  and  z  variables 

are  decoupled  the  off-diagonal  elements  of  tfie  matrix  Alt)  are 

zero.  The  width  of  the  packet  (i.e.  the  coherence  length)  in  the 

1/2 

2-direction  is  (h/ImA^^)  .  By  using  the  packets  to  construct 

the  initial  wave  function  at  t=t  we  fix  r  and  k  .  but 

leave  A_(t^).  i  =  x.y  or  z.  unspecified.  e.xcept°for  the'^  obvious 
requirement  that  if  we  fix  the  number  of  packets  used  for  the 
fit.  their  width  in  the  x  and  y  directions  should  be  large  enough 
that  they  overlap.  At  first  sight  this  feature  is  unpleasant, 
since  we  would  prefer  that  the  initial  wave  function-which 
carries  all  the  information  regarding  how  the  beam  was  prepared- 
should  fix  the  initial  conditions  for  all  the  parameters  in  all 
the  Gaussians.  Since  these  parameters  are  propagated  in  time  by 
first  order  differential  equations,  their  final  values  depend  on 
the  initial  ones.  There  is  a  danger  that  by  making  use  of  the 
existing  flexibility  in  the  choice  of  A..(t  )  we  can  affect  at 

1  J  0 

will  the  final  wave  function  (i.e.  the  diffraction  results). 

1  2 

However,  Heller  '  has  made  an  ingenious  use  of  this  freedom  to 
optimize  the  calculation:  he  choses  A. .it  )  such  that  a  freely 


will  have  a  minimu.n 


propagating  packet  starting  at  r 
uncertainty  (i.e.  ReA..(t)=0)  at  the  surface  and  a  very  small 

width.  Thus  the  packet  should  be  narrow  in  the  region  where  the 
potential  varies  fastest,  satisfying  the  dynamic  requirement  of 
the  SHM  at  the  point  where  they  matter  .most.  The  flexibility  is 
thus  greatly  diminished  but  not  eliminated.  One  wonders 
therefore  whether  the  choice  for  ImA^^  affects  the  final  results. 
Furthermore,  it  is  not  clear  that  by  choosing  the  initial 
conditions  in  this  way  we  can  control  the  width  in  the  collision 
region.  These  questions  are  answered  in  what  follows  by  carrying 
out  numerical  calculations. 


We  consider  the  case  of  He  scattered  by  a  potential  that 
mimics,  very  crudely.  Pt(lll).  The  interaction  with  the  surface 
is  taken  to  consist  of  two  parts  and  is  written  as 


V(  r  ) 
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where  R.  denotes  the  position  of  the  i-th  surface  atom.  T.he 
first  term  is  a  sum  of  Lennard-Jones  interactions  of  He  with  a 
small  group  of  surface  atoms.  This  term  is  responsible  for  the 
x-y  variation  of  the  potential  and  was  chosen  to  give  a 
corrugation  of  .19  A  for  a  He  beam  energy  of  1.2  kcal  mole  The 

second  is  a  bulk  term  which  results  from  the  pairwise  summat  .on 
of  the  Lennard-Jones  interaction  with  all  of  the  atoms  in  the 

_  3 

solid.  It  is  only  a  function  of  z.  has  the  proper  z  long  range 
attractive  part,  and  was  made  to  produce  a  well  depth  of  12. a 
kcal.  mole.  This  particular  form  for  the  potential  was  chosen  not 
to  mimic  the  actual  He,'Pt(lll)  system,  which  has  negligible  off 
specular  s ca t t e r i n g , ^ ^  but  to  provide  a  model  system  with  a 
reasonable  corrugation. 


In  Fig.  1  we  plot  the  width  of  a  packet  in  the  z  direction, 
j  as  a  function  of  the  distance  from  the  surface  for  various 


initial  values  of  A(o)  chosen  such  that,  if  we  assume  a  free 
space  propagation,  the  packet  would  be  narrowest  at  the  turning 
point.  This  can  be  achieved  if  we  determine  A(o)  in  the 


following  manner.  Since  in  free  space  there  is  no  coupling 
between  the  x . y .  and  z  coordinates,  A(t)  remains  diagonal,  and  we 
can  simply  consider  the  propagation  of  some  diagonal  component  A 
given  by 


Alt) 


2  i 
m 


i  A 


T _ 

(  t-T) 


1 


At  the  turning  point.  t=T,  A(T)  s  iA.^,  and  the  packet  reaches  its 

minimum  width  at  the  collision  site.  Thus  A(o)  is  determined  by 

using  the  time  T  and  the  chosen  value  for  A.j.  in  the  above 

equation.  In  Fig.  1  we  make  all  3  diagonal  elements  in  Alol 

equal,  and  choose  A.p  =  .01,  1,  and  100.  corresponding  to  free 

sp.ace  Gaussian  widths,  at  the  classical  turning  point,  of  2.5. 

-  1  4 

25,  and  .025  A  respectively.  The  units  of  a.^.  are  amu  (10 

sec).  The  potential  is  plotted  directly  below  in  Fig.  2  for  easy 

comparison.  We  observe  a  free  space  propagation  up  to  2  A  or  so 
2  2 

when  the  3  V.'3z  term  in  the  SHM  equations  starts  modifying  the 

width  and  gives  it  values  differing  from  the  free  space  ones.  In 

all  cases  the  broadening  is  out  of  control;  note  that  in  all  the 

trajectories  the  dynamic  effects,  present  (in  SHM)  through 
2  2 

3  V  3z  .  dominate  the  time  evolution  of  the  width  and  broaden  it 
to  about  the  same  width,  regardless  of  our  choice  for  a.^.  ;  this 
width  e.xceeds  the  values  permitted-  by  the  validity  condition  of 
SHM.  Thus  due  to  these  dynamic  effects  the  interaction  with  the 
lattice  increases  the  beam  coherence;  unfortunately  the  increase 
is  so  large  that  the  validity  of  SHM.  used  in  the  propagation,  is 
uncertain . 

In  Table  1  we  list  the  peak  heights  for  the  specular  and  a 
few  low  order  diffraction  peaks,  calculated  using  various  initial 
values  for  Alt)  as  above.  We  choose  A-  =  .1,  1.  10.  100.  1000. 


and  normalize  such  that  the  height  of  the  specular  peak  (0.0)  is 
1.0.  As  mentioned  previously  we  see  a  non - n eg  1 i g  i  b 1 e  dependence 
on  the  parameter  A.^, .  The  results  converge  for  A.^,  >  50  or  so.  In 
this  regime  a  is  very  broad  over  most  of  the  trajectory,  making  a 
quadratic  expansion  of  the  potential  invalid.  The  above  analysis 
has  also  been  applied  to  a  similar  system  with  a  purely  repulsive 
exponential  gas  surface  interaction  of  the  form 

_  -  -B  "r -R . 

V  (  r  )  =  Z  C  e  1 

i 

where  R.  again  labels  the  position  of  the  surface  atom.  We  chose 
^  -  1 

B  =  2.1  A  and  selected  c  to  give  a  corrugation  similar  to  that 
of  the  previous  example.  This  interaction  is  softer  and  of 
longer  range  than  the  Lennard- Jones  potential  used  above  but 
leads  to  similar  results  and  conclusions. 
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IV.  DISCUSSION’ 

The  analysis  of  the  preceding  section  used  the  si. n pie  casf* 
of  neutron  scattering  to  establish  a  minimal  set  of  requirements 
for  the  use  of  wave  packets  to  construct  diffracted  wave 
functions.  As  expected  the  requirements  on  the  dependence  of  "he 
incident  wave  function  on  the  "parallel"  variable  .x  and  y  are 
quite  demanding.  For  the  calculation  of  scattering  by  ordered 
systems  at  T=0  K.  or  for  the  calculation  of  the  diffraction  peak 
intensity  decrease  in  elastic  scattering  from  ordered  lattices  at 
finite  temperature,  or  for  scattering  from  disordered  lattices 
lacking  spatial  correlations,  the  "parallel"  coherence  area  of 
the  beam  must  exceed  the  area  of  the  unit  lattice.  Bloch's 
theorem  can  then  be  used  to  construct  the  scattered  wave  function 
for  the  whole  lattice.  Inelastic  scattering  from  ordered 
lattices  at  finite  temperature  or  from  disordered  lattices  with 
spatial  correlation  require  a  parallel  be. am  coherence  area  equal 
or  larger  than  the  coherence  area  of  the  density -  density 
correlation  function  of  tlie  two  dimensional  scattering  centers 
These  demands  reflect  the  fact  that  if  the  beam  is  to  be  used  to 
study  the  surface,  its  coherence  must  exceed  "he  surface 


conerence 


If  this  is  not  the  case  we  end  up  using  "he  surfitr^ 


to  study  the  coherence  properties  of  t!ie  bea.m.  u  t'lis  :s  .in- 
one  of  the  aims  of  surface  science. 


The  only  existing  comparison*  between  an  S’riM  diffraction 
calculation  and  the  corresponding  "exact"  quantum  results 
i obtained  by  solving  the  Schroedinger  equation  numerically  by 
isir.g  the  standard  eigenfunction  expansion  method)  shows  good 
agreement.  The  numerical  experience  (the  use  of  the  experiment.!  i 
uata  to  test  the  reliahiiity  of  the  method  is  not  possible  since 
t)o  little  is  known  about  the  potentials)  accumulated  so  fir  is 
■>')  limited  "o  allow  us  to  say  whether  ti\is  is  an  accident  or 
V*  -1  “  t  h  e  r  diffraction  is  peculiarly  insensitive  to  the  errors  mad- 
oy  SHM.  From  a  practical  point  of  view  we  cannot  be  satisfied 


.*• 


with  a  new  method  whose  reliability  can  only  be  estaolished 
through  a  comparison  with  the  results  obtained  wit  I',  the  me  thou 
which  it  intends  to  replace;  furthermore,  for  the  most  i. m  porta  r.  : 
applications,  to  disordered  surfaces,  we  cannot  benefit  from  sjt 
a  comparison  since  the  traditional  methods  are  ;,m  practical  for 
this  case.  It  is  therefore  of  importance  e;\umi;ie  the  fears 
raised  by  the  use  of  SH.M  and  the  improvements  that  will  aiieviat 
t  h  o  m  . 

'  1  )  The  choice  of  the  initial  state.  We  find  that  i a  )  t  h 
use  of  a  narrow  wave  packet  for  the  z-dependence  of  ti'.e  wave 
function  is  very  likely  to  affect  diffraction  intensity.  To 
prevent  this  the  packet  must  be  broader  than  the  lattice  spacing 
de  Broglie  wave  length  of  the  atom  and  the  range  of  the 
potential.  This  suggests  that  we  should  work  wi^h  broad  packets 
and  therefore  abandon  the  expansion  of  the  potential  lused  in 
SH.M)  around  the  center  of  the  packet  which  is  valid  only  wnen  th 
packet  stays  narrow  throughout  the  collision.  ibl  Another 
troublesome  feature  is  the  choice  of  the  width,  to  give  a  narrow 
mini, mum  uncertainty  wave  packet  at  the  wall,  which  is  not  unique 
does  not  have  the  effect  intended  (since  the  interaction  broaden 
•lie  packet  anyway)  and  it  is  not  physically  j  u  s  t  .  :  :  a  u  1  e  .  bur 

I  ! 

application  of  SHM  to  curve  crossing  *  or  to  the  propagition  if 
1  2 

Morse  states  reaches  the  same  conclusion.  There! a  "e  wp  sugges 

the  use  of  an  initial  state  i{i.(x.y.z)  which  is  ?.  two  dimensi  onal 

planar  wave  m  u  1  ,t  i  p  1  i  e  <1  t  >3  a  broad  packet  in  t e  z  -  d  1  r  e  c  t  1 13  n 

This  represents  correctly  t!ie  coiierence  of  the  beam  used  In 

e  .x  pe  r  i  m  e  n  t  a  1  studies.  We  can  then  fit  tills  state  ’with  a  sum  of 

Gauss  ians  E  G  (.x.y.z;(X  }  .  'where  (Xa  }  is  a  sy-mool  for  -he  se-  o 
■j.  a  o 

parameters  in  the  Gaussian  a.  by  minimizing  'D.  -  E  w;:h 

resoect  to  all  X  Thus  the  initial  values  of  a .  1  r  ,i  r  amp  t  e  :■  s  'i  r 
■  a 

fined  by  the  physical  conditions  of  the  problem.  We  h  a e  found 
tr.  .at  this  procedure  improves  d  rama  t  i  c  a  1  i  y  ‘  "  the  results  in 
;;ii'..pr  prohlp-n  that  Is  troublesome  for  SilM  .  ‘hat  d:  the 
r  '.  p  a  g  a  t  I  3  n  of  Morse  eigenfunctions.*"’ 


( 2  )  The  propagation  method.  S H M  can  be  replaced  by  a 
vMriational  method'^  which  generates  coupled,  non-linear  first 
order  differential  equations  for  Che  parameters  Thus  the 

parameters  in  the  Gaussian  a  are  affected  by  the  motion  of  all 
other  Gaussians  p  through  coupling  terms  that  become  non-zero 
whenever  at  and  3  overlap.  The  presence  of  this  coupling  is 
clearly  a  necessary  and  important  featufe  since  if  two  Gaussians 
are  present  in  the  same  spatial  region  SiT)  centered  around  T 
they  construct  the  scattered  wave  function  i  r  ;  t  )  for  reS(r) 
together  :  it  is  inconceivable  that  they  can  do  this  with  any 
accuracy,  if  they  propagate  at  all  time  independently  of  each 
other.  This  "inCergaussian"  coupling  affects  all  the  parameters 
■  iie  width,  the  phase,  the  mean  position  and  the  mean  momentum  of 
the  packets.  It  changes  the  physical  picture  from  that  of  a  set 
of  indepenent  Gaussians  moving  along  classical  trajectories  to  a 
‘community"  of  Gaussians  moving  on  non-classical  trajectories 
that  are  coupled  to  eacli  other. 


Furthermore  there  is  no  need  and  no  general  pliysicol  or 
mathematical  justification  for  expanding  the  potential  around  th 
center  of  the  Gaussians.  If  the  potential  is  fitted  to 
exponentials.  Gaussians  or  polynomials,  or  any  combination  of  'h 
aijove.  all  the  integrals  appearing  in  the  theory  can  be  done 
analytically  and  very  little  of  the  SHM  simplicity  is  lost.  One 
of  the  interesting  consequences  of  giving  up  the  e.xpansion  is 
that  the  force  3V(r  )  3r  appearing  in  SHM  is  replaced  by  tiie 
derivative  of  W(r^l  =  fG  (r)V(r)G{r)dr  and  this  time  dependent 
potential  differs  in  interesting  ways  from  V(r^)  For  example'^ 
as  the  packet  approaches  the  surface  W(r^)  becomes  non-zero  .as 
Soon  as  the  Gaussian  overlaps  with  the  potential 

differs  from  the  classical  potential  V(r.i.  claanges  in  time,  and 


'-‘nds  to  come  forward  and  meet  the  packet.  The  packet  turns 
iriumi  earlier  than  in  the  classical  case  (i  e.  SHMl.  Since  the 
mo  "ion  of  “he  center  of  the  uncoupled  packets  takes  place  on 
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•»(r^:t)  it  is  conceivable  that  this  has  a  different  corrugarion 
th(in  the  true  potential  and  thus  affects  diffraction. 

We  are  currently  pursuing  these  improvements  of  SHM  For 

11  ■  ^ 

Che  case  of  curve  crossing  and  of  the  .Morse  oscillato:"  they 

turn  out  to  be  essential.  For  diffraction  we  are  .not  yet  sure. 

1  2 

The  striking  success  of  SHM  '  is  partly  due  to  Heller's 
e.\c  client  choice  of  problems,  which  involve  very  iocaltzed  wave 
functions  that  need  to  be  propagated  for  a  very  short  time  anc  or 
harmonic  oscillator  problems  in  which  the  Gaussiar. s  do  not 
o’‘erlap.  Diffraction  deals  with  a  very  delocalised  wave  function 
in  which  the  overlap  between  different  Gaussians  is  very  likely. 

Thus  it  is  difficult  to  untlerstand  why  SHM  should  work  well  ;n 
this  case.^  A  conceivable  reason  might  be  that  the  process  is 
dominated  by  interference,  which  SHM  treats  at  a  semi  - classical 
level  (since  it  gives  each  trajectory  a  phase  equal  to  the 
classical  action). 

For  problems  in  which  the  simplifying  features  .mentioned 
above  (i.e.  localization  and  short  time  processes)  are  absent, 
more  laborious  versions  of  SHM  might  be  needed.  Our  limited 
experience  indicates  that  even  with  tlie  computational  effort 
added  by  these  improvements  Heller’s  method  still  has  e.xtremeiy 
good  prospects  to  contribute  substantially  to  surface  science 
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